Dirac-like approach for consistent discretizations of classical constrained theories 



Cayetano Di Bartolo 1 , Rodolfo Gambini 2 , Rafael A. Porto 3 and Jorge Pullin 4 
1. Departamento de Fisica, Universidad Simon Bolivar, 
Aptdo. 89000, Caracas 1080-A, Venezuela. 
2. Instituto de Fisica, Facultad de Ciencias, 
Igud 4225, esq. Mataojo, Montevideo, Uruguay. 

3. Department of Physics, 
Carnegie Mellon University, Pittsburgh, PA 15213 
4- Department of Physics and Astronomy, 
Louisiana State University, 
Baton Rouge, LA 70803-4001 

We analyze the canonical treatment of classical constrained mechanical systems formulated with 
a discrete time. We prove that under very general conditions, it is possible to introduce nonsingular 
canonical transformations that preserve the constraint surface and the Poisson or Dirac bracket 
structure. The conditions for the preservation of the constraints are more stringent than in the 
, continuous case and as a consequence some of the continuum constraints become second class upon 

discretization and need to be solved by fixing their associated Lagrange multipliers. The gauge 
invariance of the discrete theory is encoded in a set of arbitrary functions that appear in the gener- 
ating function of the evolution equations. The resulting scheme is general enough to accommodate 
the treatment of field theories on the lattice. This paper attempts to clarify and put on sounder 
footing a discretization technique that has already been used to treat a variety of systems, including 
Yang-Mills theories, BF-theory and general relativity on the lattice. 
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I. INTRODUCTION 



We have recently introduced a technique for treating the theories that arise when one discretizes (space)- 

time in a constrained mechanical system or a continuum field theory. We have shown that this technique works for 
Yang-Mills and BF theories and implemented it for the gravitational case. Previous attempts to studying systems 
with discrete time had concentrated mostly on systems without constraints or with holonomic constraints only (for a 
| review with a comprehensive reference list see Q.) 

The idea consists on starting from a discretized action, constructing discrete Lagrange equations and introducing 
a symplectic structure in the discrete space. The evolution is implemented via canonical transformations and the 
consistency of the discrete theory determines in part the Lagrange multipliers. In some totally constrained systems, 
^Jy like general relativity, the resulting discrete theories are constraint-free since the constraints are solved for the Lagrange 
multipliers. This makes the quantization of the discrete theories considerably simpler than the continuum cases. This 
was exploited to make progress in solving the problem of time in quantum gravity 0, and to implement the Page- 
Wootters relational time and show that a fundamental decoherence arises in quantum mechanics from quantum 
gravity. 

In this paper, we want to address in a more systematic way the issue of the canonical formulation of discrete 
constrained systems. Up to now, most of the analysis has been made on specific examples, and a canonical analysis, 
a la Dirac0,13, ^ s still lacking. In particular, the technique relied heavily on defining a canonical transformation 
that was initially singular (and therefore not a true canonical transformation) and showing that one could eliminate 
variables and end up with a true canonical transformation. Up to now this was shown in a case by case basis. This 
paper determines the general conditions needed for the construction of a proper canonical transformation by following 
a close analogue of the Dirac procedure, adapted to the discrete case. In particular, we note that there are several 
ways to proceed that yield equivalent results, but that may offer different advantages for particular systems. 

In section II we will lay out the framework of how to deal with mechanical systems where time is discrete, including 
singular and non-singular systems. In section III we will develop a classification of constraints into first and second 
class suitable for the discrete context. In section IV we work out a specific example that exhibits the details of the 
formalism. In section V we discuss an alternative formulation of the formalism and we end with conclusions and 
discussion. 



II. MECHANICS WITH DISCRETE TIME 



We start by considering a continuum theory representing a mechanical system. Its Lagrangian will be denoted by 
L(q a , q a ), a = 1 . . . M. This setting is general enough to accommodate, for instance, totally constrained systems. In 
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such case q will be the derivative of the canonical variables with respect to the evolution parameter. It is also general 
enough to include the systems that result from formulating on a discrete space-time lattice a continuum field theory. 

We discretize the evolution parameter in intervals (possibly varying upon evolution) t n +i — t n = e n and we label 
the generalized coordinates evaluated at t n as q n . We define the discretized Lagrangian as 



L(n, n + 1) = L(<£, q a n+1 ) = e n L(q\ q a ) 



(1) 



where 



q = q n and q = 



(2) 



The action can then be written as 



S 



N 

£ 

n=l 



L(n, n + 1) 



(3) 



If the continuum theory is invariant under reparameterizations of the evolution parameter, one can show that the 
information about the intervals e n may be absorbed in the Lagrange multipliers. In the case of standard mechanical 
systems it is simpler to use an invariant interval e n = e. 

The Lagrange equations of motion are obtained by requiring the action to be stationary under variations of the 
configuration variables q a fixed at the endpoints of the evolution interval n = 0, n = N + 1, 



dL(n,n + l) dL(n — l,n) 



= 0. 



(4) 



These equations define a unique evolution if the determinant 

d 2 L(n,n+ 1) 



(5) 



We will refer to this case as the nonsingular case. When the determinant vanishes, one will have to analyze the 
situation differently. Let us start with the nonsingular case. 



A. Nonsingular case 

In this case one can solve the Lagrange equations explicitly and the q n +i are uniquely given as a function of q n and 
q n -i- This is the equivalent of the Hessian condition for the non-singular Lagrangian theories in the continuum. The 
resulting equations are "second order" in the sense that the g„'s are determined provided one knows two previous 
time levels. One can introduce a "first order" formulation by introducing canonically conjugate variables as is usually 
done when introducing a Hamiltonian formulation in the continuum theories. 

We introduce the following definition of variables that we will later show end up being canonically conjugate 
momenta of the configuration variables, 

dL(n,n+l) 

Pn+l ~ g a W 

a = dL(n - l,n) ^ dL(n,n+ 1) 
Pn ~ dq° dq- [ > 

Where we have used Eq. (0J. The equations (JSJ and Q define a canonical transformation for the variables q n ,Pn to 
q n +i,Pn+i with a the type 1 generating function Fi = —L(q!^,q^ + 1) provided that condition Q is fulfilled. Notice 
that the evolution scheme is implicit, one can use the bottom equation (since we are in the non-singular case) to give 
an expression for q n +\ in terms of q-mPm which in turn can be substituted in the top equation to get an equation for 
p n+ i purely in terms of q n ,p n . 

It should be noted that there are several other possible choices, when going from the set of equations I|6I7|I to an 
explicit evolution scheme. For example, one can choose to do things in a way that yields a closer analogy with the 
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standard Hamiltonian description in the continuum by introducing type two canonical transformations. To do this, 
we choose to invert equation © for q% +1 = qn+i(ln'Pn+i)> which is possible only if 



d 2 L{n,n+ 1) 



We can now introduce a Legendre transform and define 



^0. 



F 2 {q n ,p n+1 ) = ^CflPn+l ~ L (5n,en+l) 



(8) 



(9) 



From here it is immediate to obtain 



Hn+1 



dF 2 
dq b n 



dF 2 
dp b n+ i 

-dL(n,n + l) 
dql 



Pi 



(10) 

(11) 



where we have used the Lagrange equation in order to obtain the last equality. We easily recognize from here that 
F 2 behaves as a type two generating function of the canonical transformation connecting level n with level n + 1. 

We may define now a sort of "type 2" Hamiltonian (in the sense that it depends on p n +i and q n , H 2 (j> n+ i, q n ) given 
by 



H 2 (q n ,p n+1 ) = 2jK +1 (g« +1 ~ In) ~ L (ln,q n +i) = F 2 (q n ,p n+1 ) 



.Pn+iq-n 



(12) 



which leads to the discrete Hamilton equations, 



Hn+1 



Pn 



dH 2 



d P b 



Pn 



71+1 

dH 2 
dq b n ' 



(13) 
(14) 



It should be noted that although this formulation has a degree of analogy with the traditional Hamiltonian for- 
mulation, there are significant differences due to the fact that the conjugate variables live at different time slices. It 
would not be possible therefore to use this formulation to attempt to construct a Schrodinger equation starting from 
the above Hamiltonian. 

Provided that the canonical map defined by F 2 is invertible we end up with a discrete evolution implemented by a 
canonical transformation. It can be easily seen by using the Legendre transform that F 2 will be invertible if and only 
if 



d 2 F 2 (n,n + l) 



dq*p b 



nfn+l 



d 2 L(n,n+l) 



dqvqf 



n+l 



d 2 L(n,n+l) 



-1% 



?0. 



(15) 



Thus, in order to have canonical transformations generated by type two functions, the Hessian condition in the 
continuum time mechanics leads to two independent conditions in the discrete theory given by Eqs. © and (JHJl. 
Notice however that JSJ) is not necessary for introducing a symplectic structure. 

It is clear that when one builds a canonical discrete theory there are four possibilities depending on which pair of 
variables one chooses to construct the generating functional of the canonical transformation, either q ni q n +i, 9n,Prt+i, 
p n ,q n+ i, p n ,p n+ i. In this subsection we considered only two, but the others can be easily generalized from the 
discussion here. 



B. The singular case 

Let us consider as before the generic discrete Lagrangian L(n, n+l) = L(qf l7 q„+i) with a = 1 . . . M. It leads to 
the equations we already discussed, 



dL(n, n+l) 
dqJTx 



PUi = ' (16) 



Pn - -^T^ d7) 
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If 



d 2 L{n,n+ 1) 



vanishes and its rank is K, the system is singular and has M — K constraints of the form 

<MCK) = o (is) 



that result from Ea (|17|l . and M — K constraints 



MC+nK+i) = (19) 

resulting from Ea H16() . The evolution of the configuration variables from level n to n + 1 is given by solving for 
the Eqs. p7jl. As the system is singular the evolution depends of (M — K) arbitrary functions V A 

q a n+1 = t{qtplV A ). (20) 

We shall follow closely the standard Dirac canonical procedure of continuum mechanics. In this case, the analysis of 
a constraint system goes trough two steps. The first step consists in the definition of a set of evolution equations that 
weakly preserve the constraints and the Poisson symplectic structure. To do that one defines the total Hamiltonian 
Hf = Hq + V a <j) a where <p a are the primary constraints, and the V a are partially determined in order to preserve all 
the constraints of the system. Even though some of the V's may be arbitrary functions, once they are specified, the 
evolution generated by Ht preserves the Poisson brackets and the dynamical evolution is consistent at the classical 
level with the constraint structure. The second step is only required to quantize the system and consists in the 
identification of the first and second class constraints and the introduction of the Dirac brackets that enforce strongly 
the second class constraints. 

As we shall see, the same procedure may be followed in the discrete case. The main difference is the implementation 
of the canonical transformation that is not generated by a Hamiltonian but by a canonical transformation of type 2, 
3, or 4. 

Let us start by completing the evolution equations. We need to add to Ea (|20|) an equation for p n +i 



a _ dL{q n ,q n+1 ) 

Pn+l — 



dQn+1 

We now impose the preservation of the constraints 



(21) 



W n+1 ,p a n+1 ) = ^A ^ J^p ) =0. (22) 

Furthermore, we need to impose the n + 1 level constraints at level n, ^A(qn,Pn) = and impose the consistency 
conditions 

Three different cases may occur: 

a) Eas. l|22l23|l vanish automatically, and therefore we are not led to new conditions. 

b) They lead to inconsistencies, and the dynamical system is inconsistent. 

c) New secondary constraints C(q^,p^ l ) appear or/and some of the arbitrary functions V A are determined, that is 
V A — V A (q^,pf l ,v a ) with a = 1...R < (N — K), and v a arbitrary functions. The process is repeated until 
consistency is achieved. That is, until the consistency conditions are automatically satisfied without further 
constraints and conditions for V. 

Substituting V A in l|20|l and l|21|l we get the evolution equations that preserve all the constraints: primary, secondary, 
tertiary, and so on, 

= f a (ql P tV A (q,p,v)) = f a (ql P tv a ) (24) 




P a n+1 =9 a (q b n ,p b n ,v a ) 



(25) 
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Initial values need to be restricted by the n = level constraints, 

MM) = *A(<Z a ,Po a ) = 0. 



(26) 



In order to have a complete analogy with the continuum case, we still need to analyze under what conditions 
this evolution also preserves the Poisson bracket structure. As in the continuum case we assume that the arbitrary 
functions v have been fixed. Three different cases may arise depending on if one chooses to implement things in terms 
of a canonical transformation of type 2, 3 or 4: 



Case I) The equation 1(33} is invertible for q^: that is, 



^ and therefore one can write 



Qn = h a (q b n+vP h n ). 

Notice that, under these hypotheses, there are no pseudo-constraints of the form G(g^ +1 ,p°) 
pseudo-constraints because they involve variables at different instants of time. 
We may define a type 3 generating functions of canonical transformations, 



Fs(q b n+ i,p b n) 



[p\q b n + L(q n ,q n+1 )} Lb h t. 



Then we have 



and 



a*3(g£+i,p£) 



dF 3 (q? l+1 ,p«) _ dL(q a n ,q- +1 ) 



9q b n+1 



= P b n, 



(27) 

0. We call these 
(28) 

(29) 
(30) 



Qn—h 



Notice that in the last equality there are also contributions coming from the dependence on q n +i of the level n 
variables q n , but these contributions cancel because of the definition of the canonical momenta. The information about 
the momenta is completely encoded in the evolution equations q^+i = h a and the constraints. As the first equation is 
equivalent to l|24|) one ends up recovering the fundamental evolution equations as a canonical transformation generated 
by F3. Furthermore, 



&F3 Pi) 




dq b 


dp b n dql +l 




dq a n+1 



^0 



(31) 



due to the fact that we have assumed that 



dq b n +i 



and therefore the resulting canonical transformation preserves the Poisson bracket structure 

dp a 

Case II) The equation is invertible for p ^. that is "" 



^ 0, and consequently F3 is a non singular generating function 
5 Poisson bracket structure. 
7^ and therefore one can write, 



dp b n 



P a n = g a (qlp b n+1 ) 



(32) 



Notice that, under these hypotheses, there are no pseudo-constraints of the form G(q^,p^ l+1 ) = 0. By substituting 
(|32l in one gets 



which allows to introduce a type 2 generating function 



(33) 



(34) 



One can now easily check that this generating function reproduces the evolution equations (|24|l and l|25|) and defines 
a non singular canonical transformation that preserves the Poisson Brackets in the evolution. 

Case III) Even when the system has pseudo-constraints of the form G(q® , i,p%) — and F(q') ll p^ l+l ) = one may 
be able to find a canonical transformation provided that the system does not have pseudo-constraints of the form 
F(p a n , P a n+1 ) = 
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In fact, by using Ea. (|25|l one can invert for 



•n 



a 



(35) 



and substituting in ipijl get 




m a (p b n ,p b n+1 ). 



(36) 



A generating function of type 4 that does the same job that the two previous ones may now be introduced, 



Fi(p b n ,p h n+l ) = [-P b n+1 q b n+1 +P b n q b n + L(q n ,q n+1 )] 



(37) 



All of the discrete systems that have been treated up to now in the literature may be analyzed by following this 
canonical procedure, allowing to preserve the constraints and the Poisson bracket structure. Later on we will show 
an example of a system of this type in order to analyze how this procedure works in a concrete case. 

It should be noted that there may exist mechanical systems that do not fall into any of the above classifications. 
For instance, a system could have pseudo-constraints of all the types listed above. In such cases one will need to 
develop further techniques to treat them. For instance one could introduce canonical transformations of a given type 
for some of the variables and of a different type for other variables. This would require further study and it does not 
appear necessary for the systems that have been analyzed up to present. 



At this point we have a set of constraints primary, secondary, tertiary, etc, of the form Xz(q a ,p a ) — with 
Z = I... A with A the total number of independent constraints, that are preserved under the evolution given by the 
equations Ij24l25|) provided part or all the arbitrary functions V are conveniently fixed. 

As in the continuum case it is convenient to introduce the notion of first and second class constraints, in order to 
quantize the theory. A constraint is of first class if it commutes with all the constraints, if that is not the case it is of 
second class. As in the continuum case one can define first class functions of the canonical variables f(q,p) that are 
not necessarily constraints. Such a function will be first class if it commutes with all the constraints. Second class 
constraints may be imposed strongly by introducing Dirac brackets. As the evolution equations preserve the Poisson 
structure, they will preserve de Dirac structure because Dirac brackets arc defined in terms of Poisson brackets. One 
ends up with a theory with a set of evolution equations that preserve the symplectic structure of the system, and 
therefore may be quantized by describing the evolution in terms of unitary operators. 

In the discrete case there is not a straightforward relation between the number of first and second class constraints 
and the number of phase space degrees of freedom. This is due to the fact that now the evolution of the constraints 
is not directly related with their Poisson brackets with a total Hamiltonian. Thus, the fact that a constraint does 
not commute with others is not easily related with the determination of an arbitrary function. It is still very easy to 
determine the number of phase space degrees of freedom. In fact this number is given by two times the number of 
configuration variables minus the total number of constraints minus the number of arbitrary functions v. 

IV. AN EXAMPLE OF CONSTRAINED SYSTEM WITH SECOND CLASS CONSTRAINTS 

To illustrate the techniques elaborated above, we would like to discuss a model that is simple, yet addresses in 
a non-trivial way the main points we discussed. This example had been treated using ad- hoc techniques in Q. 
The model consists of a parameterized free particle in a two dimensional space-time under the influence of a linear 
potential. The discrete Lagrangian is given by, 



We have chosen a first order formulation for the particle. However, this Lagrangian is of the type we considered in this 
paper, one simply needs to consider all variables, q a ,ir a ,N as configuration variables. The system is clearly singular 
since the it's and N only appear at level n (or in the continuum Lagrangian, their time derivatives are absent). When 



III. CLASSIFICATION OF THE CONSTRAINTS 



L n = L{q a n , <, N n , < +1 ,iV n+1 ) = K(ln+i - In) - N n K + -K 1 ,) 2 + aql]. 



(38) 



7 



considered as a Type I generating function, the above Lagrangian leads to the equations 



Pln+l = ^=0, (39) 

un n+l 



dL n 



Pln+l = ^-=C (4°) 
U! ln+1 



^ = 8^ = °' (41) 



and 



Pin = " flZJ = - 9n) + + (42) 

Pin = ~=K+^ n , (43) 
<9L 1 

PN,n = -d]^^ 7T n+2^ 2+aq ^ (44) 

One can easily recognize that the system has six constraints: three at the n+ 1 level, and three at the n level. They 
are: 

ri = Pln+1=0, (45) 
V>2 = P N , n+1 = 0, (46) 

*i = p» n -« + 6?aJV„), (47) 
*2 = Pjv, n - [tt° + i(7ri) 2 + a<£]. (48) 

Therefore the evolution depends on three arbitrary functions V N n ,V® n , 

Qn+i = Qn + <N n 6 a x + N n 5 a - p« „ (49) 

< +1 = K + v?, n (50) 

iVn+i = N n +V Njn . (51) 

The preservation of the ip constraints from level n to level n + 1 is automatically ensured from (|23|) . Now we impose 
the preservation of the $ constraints upon evolution. Let us begin with 

= Pi n+1 - <+l =Pin~<~ K n = 0, (52) 

which taking into account the constraint hnplies V° n = 0. 
For the $2 one gets the equation 

<S>2n+i = tt° + a(<£ + ttX) + « + C') 2 /2 - 0, (53) 
that taking into account the constraint $2 n implies that 

Vl „ = -4 + eVWJ 2 - 2n n N n a (54) 

where e = ±1. 
Finally we have 

<*> L+i = Pi, n - N n a - n n - Vl „ - a(N n + V N> J = 0, (55) 
that after imposing the constraint at level n leads to 

V N , n = ~Vl n -N n . (56) 
a 
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Thus, the evolution equations for the configuration variables are 



Qn+i = q^ + ^ n N n 5t+N n 5^ 



'n+l 



7T„ 



4+1 = ^i<) 2 - 2<N n a 
N n+1 ' 



(57) 
(58) 
(59) 

(60) 



We are now ready to define an invertible canonical transformation with the help of a type 3 generating function. 
Notice that these evolution equations are invertible for and therefore we are in the case I. The inverse is given by 
a set of equations of the form q n (q n+ i), explicitly given by, 



( lu 



n^ +1 N n+1 + |(AWi) 2 

Hn+1 



'n+l 



aN n+ i 



'n+li 



'n+l 



aN, 



n+l, 



^n+l^n+l + TT(^»+l) 2 



"n+l f ®N n+X 

Recalling that F3 is given by equation Ij28(l we obtain, 

1 



(61) 

(62) 

(63) 
(64) 

(65) 



^3 = P°,„7T°+i +p^ n (aN n+ i + ttI +1 ) + p°, n q° +1 + T; N n+i [ot 2 Nl +l + ?>aN n+1 Tr] l+l + 2(tt^ +1 ) 2 ] 



N n+1 (aN n+1 + 2wj +1 ) 
2(aN n+1 +7ri +1 ) 



+ P\,n [Qn+l - ^ N n+1 ~ iV «+l 7r n+l) ■ ( 66 ) 



One can check that Eqs. I|61|) to (|65|l are easily recovered by taking the partial derivative with respect to By 
differentiating with respect to q^ +1 one gets 

(67) 
(68) 
(69) 



Pq, n+l 


= P° 




Pq, n+l 


= P q ,n ~ ^A n N n+1 (aN n+ i + 2tt* +1 ), 




rP 






Pit, n+l 


= P° 




P-rt, n+l 


= <„ + ^ [3«iV n+1 - 2p\ >n + 4^ +1 + 


2B n A n 


PN, n+l 


3 





+aN n+1 
where we have introduced 



-P q ,n + 37r n+i ~ ^B n A 2 n (aN n+1 + 2tt^ +1 ) 



(71) 



A„ = (a^+TT^)- 1 , (72) 

B n = P N ,n-Pq,n ~ \{^l + lf ~ «9n+l- (73) 

By substituting in Eas. l|67|l to (|71(l the variables q n +i and using the level n constraints one gets Eqs. (|39() to 141|l . 
which is a canonical transformation that reproduces, on the constraint surface, the evolution equations of the discrete 
particle. 

What remains to be done is to identify the second class constraints and impose them strongly. The complete set 
of six constraints of this model ij; and $ are second class and allow to solve for 7r a and ./V and eliminate completely 
these variables and their complex conjugates P w and Pjy. 
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One can proceed in two different ways. The first alternative is to start by observing that P n and Pjv vanish strongly, 
and then to solve for ir and N in terms of the n + 1 level variables, leading to 

< = P° g ,n+l, (74) 
n n = Pq, n+lJ (75) 

N n = (76) 

"Pq, n+1 

where C„+i = p° n+1 + (p 9 ,„+i) 2 /2 + aq n+ i- The relevant evolution equations are obtained from (|57l67l68() .and are 

given by 

r< 

(77) 

(78) 
(79) 

-^ ±L , (80) 

Pg, n+1 

and we recover the evolution equations obtained in Q . 

The second alternative consists in solving for ir and N in terms of the n level variables, leading to 

< = Pin, (81) 



</n+l 


= € + 


N n = q n + 


C n +i 

UPq, n+1 ' 


Qn+l 


= Qn + 


NnKn = qi 


Cn+1 

3 

a 


P° 

fq. n 


= p° 






Pq. n 


= Pg,n- 


fl + aN n = 


Pq, n+1 + 



^-2(pO n + a q l), (82) 
= ^(p 9 ,„-ey / -2(p0„ + aq^, (83) 



and from here, computing the evolution equations by using H57I67I68JI . The two methods yield evolution schemes of 
different functional form since one propagates "forward" in time and the other "backward" . The inequivalence in the 
functional form stems from the fact that the discretization of the time derivatives chosen in the Lagrangian is not 
centered. It should be emphasized that if one starts from given initial data and propagates forwards with the first 
system of equations and then backwards using the second, one will return to the same initial data. 

Notice that we have 6 second class constraints, and the initial number of phase space degrees of freedom was 10. 
By noticing that there are no arbitrary functions left, one is left with 4 degrees of freedom on the constraint surface. 
The continuum model had two degrees of freedom. 

The procedure we have followed here is completely general and may be simplified when one is treating specific cases. 
For instance, as it happens in the continuum theory Q it is sometimes possible to implement the canonical analysis 
by first solving the constraints for the unphysical degrees of freedom N, Tr,p ,p* and then introducing a generating 
functional on the physical degrees of freedom by following the procedure of the previous sections. In this particular 



case it is easy to show that, for p q ^ 0, 



7^ 0, where q a = (q 1 ^ ), and thus it is possible to construct an 

vq n 

F3 generating functional. In all the models treated up to now in the literature the unphysical degrees of freedom 
were eliminated before obtaining the canonical transformation for the evolution of the physical degrees of freedom. 
In order to keep the analysis general in a simple model, here we have kept all the variables involved in this approach. 



V. TREATMENT IN TERMS OF TYPE II GENERATING FUNCTIONS 



Up to now, we have taken as the starting point a singular Type I generating function given by F\ = — L(q%, q^+i)- 
It is interesting to analyze how singular systems may be described in terms of other types of generating functions. 

Let us assume that neither condition JSJ nor condition (J5J arc fulfilled by the discrete system and, therefore, we are 
in the singular case. Then, it may be immediately seen, by using Ultjfl .and taking into account that the determinant 
that appears in JSJ now vanishes, that singular systems with first order Lagrangians in the continuum lead to the 
presence of pseudo-constraints 4> a {.Q.niPn+i) — 0, in the theory. Recall that these pseudo-constraints are similar to 
the constraints that arise in the continuum case, but mixing configuration variables at one level with momenta in the 
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next level. As in the continuum case, one can then introduce the pseudo-constraint surface which we shall assume 
has well defined functions, 



where a vector f in the tangent space of S r j )1 Tt q ^\S$, is such that 

t • V(f> a = Va. 
Let us introduce now as before the type-2 Hamiltonian 1 11 21 1 . 



(84) 



(85) 



(86) 



For which it is easy to see, using the pseudo-constraints, that it is a function of QmPn+i- Let us consider now an 
infinitesimal variation of H 2 along the pseudo-constraint surface, 



dHi 



E 

a 

E 



(<£+! - C)dpl+1 + [Pn+l ~ 
iln+1 - 1n)dp a n+1 - K+l + 



dL 

dq a n+1 

dL 



dL 

dq a n+1 - ( K+l + ) dq a n 



(87) 



where in the last step we used ©■ H 2 is well defined in but can be extended to the whole phase space as it is 
done in the case of constraints in the continuum theory. In order to obtain the canonical equation of motion we start 
from the identity, 



(dH dH \ 

~dr dq " 1 + dp^~ 1 dp ' l+1 ) ' 



and evaluating it for an infinitesimal displacement (dq,dp) in T^ p^S^, we can use (|87[1 and the Lagrangian equations 
of motion (0 to obtain, 



E 



Pn+l Pn + 



djh 



dq a n 



C+i + 



dH 2 
dp a n+ i 



dPn+l 



(89) 



Since (dq, dp) is an arbitrary tangent vector to S^, following (|85fl we obtain that the coefficient must be proportional 
to the gradient. Introducing therefore the Lagrange multipliers A a as the proportionality factors, we will end with 
the set of equations, 



9n+l 



ql 



dH 2 



Pn+l 
<j>a(q n ,Pn+l) = 



Pn K~„ A r, 



(90) 

(91) 
(92) 



This system is similar to the continuum set of equations with the inclusion of pseudo-constraints, which have the 
same functional form as the continuum constraints, but involve configuration variables at level n and momenta at 
level n + 1, instead of n level variables. From this initial evolution equations one may follow a procedure similar 
to the one developed in the previous section in order to study the consistency of the (primary) constraints and 
pseudo-constraints . 



VI. DISCUSSION AND CONCLUSIONS 



We have provided a canonical procedure for the introduction of a preserved symplectic structure in discrete con- 
strained systems. The analogy with Dirac's procedure in the continuum is quite remarkable. It is possible to define 
a notion of discrete evolution that weakly preserves constraints and Poisson brackets. The distinction between first 
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and second class constraints is still useful and when second class constraints are imposed strongly the resulting Dirac 
brackets are preserved. 

A feature of the discretized theories is that they may have a smaller number of first class constraints, and conse- 
quently more degrees of freedom than the continuum counterparts. The extra degrees of freedom come from the fact 
that the discrete theories may not necessarily have the same symmetries as the continuum theories. For instance, in 
the case of homogeneous cosmologies studied in 9] the extra pair of phase space degrees of freedom are associated 
with the fact that in the discrete theory different choices of refinements in the discretization in time correspond to 
different solutions in the discrete theory that nevertheless approximate the same solution in the continuum theory. 

This is only a first step for a complete understanding of the dynamics of discrete gauge systems. The relation 
between the number of constraints of first and second class and the number of degrees of freedom, and the connection 
between the first class constraints and the gauge invariance of the discrete dynamical system need to be further 
studied. Moreover, as discussed in the body of the article, if one wishes to consider more pathological systems than 
the ones considered here, more elaborate canonical transformations may need to be introduced. 

The issue of the continuum limit is well understood in the non singular case, where there is an external step 
parameter that controls the approximation. However, it needs further study in the case of singular systems, particularly 
in the case of totally constrained systems where the step of the approximation is encoded in the additional degrees 
of freedom of the discrete theory. This issue has been studied in several models0,0| but a complete characterization 
of the possible behaviors is still lacking. A similar comment applies to the role of spatial discretizations when one is 
considering lattice field theories. 
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